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(57) Abstract 

The invention relates to a faulty technical system, which is represented by an implicit stochastic differential equation system. An 
approximate solution of said system is determined by means of a discrete approximation process conforming to rule (a). The fault (b) is 
determined through iterative solving of the approximation process. 

(57) Zusammenfassung 

Es wird ein technisches System, welches einer Stoning unterliegt, mittels eines impliziten stochastischen Differentialgleichungssystems 
beschrieben. Eine naherungsweise Losung des Systems wird ermittelt, indem ein diskreter NaherungsprozeB realisiert wird. Der diskrete 
NaherungsprozeS wird gemaB der Vorschrift (a) realisiert. Durch iterative Losung des Naherungsprozesses wird die Stoning (b) ermittelt. 
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Beschreibung 

Verfahren und Vorrichtung zur Ermittlung einer Storung eines 
technischen Systems 

Die numerische Simulation elektrischer Schaltungen hat in den 
letzten Jahren grofie Bedeutung bei der Entwicklung von Compu- 
terchips erlangt. Aufgrund der hohen Kosten ftir eine Mu- 
steranfertigung eines Chips und eines moglichen Redesigns 
sind Simulatoren unabdingbar geworden. Mit ihnen ist es mog- 
lich, pradiktive Aussagen tlber Betriebsverhalten und Effizi- 
enz der modellierten Schaltung an einem Rechner zu erhalten. 
Nach erfolgreichen Simulationsergebnissen wird Ublicherweise 
ein Chip in Silizium gebrannt. 

Unter Verwendung des sogenannten Netzwerkansatzes wird eine 
Schaltung durch ihre topologischen Eigenschaf ten, die charak- 
teristischen Gleichungen von Schaltelementen und die Kirch- 
hoffschen Regeln beschrieben. 

Zur Analyse einer Schaltung wird die aus [1] bekannte soge- 
nannte modif izierte Knotenanalyse eingesetzt. Diese filhrt auf 
ein dif f erential-algebraisches Gleichungssystem der Form 



Allgemein ist unter einem dif f erential-algebraischen Glei- 
chungssystem ein Gleichungssystem des Typs 



zu verstehen mit einer singularen Jacobi-Matrix F. der par* 



C(x(t)) ■ x(t) + f (x(t)) + s(t) = 0 . 



(1) 




(2) 



x 



tiellen Ableitungen von F nach x(t) . 
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Das dif ferential-algebraische Gleichungssystem (1) wird auch 
als quasilineares-implizites Gleichungssystem bezeichnet. Mit 
x(t) wird ein von einer Zeit t abhangiger Verlauf von Knoten- 

spannungen bezeichnet, mit x(t) dessen Ableitung nach der 
5 Zeit t. Mit f(x(t)) wird eine erste vorgegebene Funktion be- 

zeichnet, welche Leitwerte und nicht-lineare Elemente ent- 
halt, mit c(x(t)) eine vorgegebene Kapazitatsmatrix und mit 
s(t) eine zweite vorgegebene Funktion, die unabhangige Span* 
nungsquellen und Stromquellen enthalt. 

10 

Die Gleichung (1) beschreibt nur den Idealfall fur eine 
Schaltung. In der Praxis lalit sich jedoch Rauschen, d.h. eine 
Storung der Schaltung, nicht vermeiden. Unter Rauschen ver- 
steht man eine ungewollte Signalstdrung, die beispielsweise 
15 durch thermische Effekte oder die diskrete Struktur der Ele- 
mentarladung verursacht wird. Aufgrund der fortschreitenden 
Integrationsdichte integrierter Schaltungen wachst die Bedeu- 
tung der prSdiktiven Analyse solcher Effekte (Rauschsimulati- 
on) . 

20 

Bei der Analyse einer Schaltung unter Berttcksichtigung von 
Rauschen kann die Gleichung (1) nodelliert werden als: 

C(x(t)) • x(t) + f(x(t)) + s(t) + B(t, x(t)) • v((D, t) = 0 . (3) 

25 

Hierbei bezeichnet v(o, t) einen m-dimensionalen Vektor, des- 
sen unabhangige Komponenten verallgemeinertes weilies Rauschen 
sind, wobei m die Zahl der Storstromquellen angibt. Eine auch 
als Intensitatsmatrix bezeichnete Matrix B(t, x(t)) hat die Di- 

30 mension n x m. Tritt nur thermisches Rauschen in linearen Wi- 
derstanden auf, so ist sie konstant. 

Wird die Schaltung mit rein linearen Element en modelliert, so 
wird die Vorschrift (3) zu einem gestorten linear-impliziten 
35 dif ferential-algebraischen Gleichungssystem der Form 
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C • x(t) + G • x(t) + s(t) + B(t, x(t)) • v(ffl, t) = 0 



(4) 



10 



15 



Unter dem Begriff Index ist im weiteren ein Mali zu verstehen, 
wie weit sich ein dif ferential-algebraisches Gleichungssystem 
von einem expliziten gewdhnlichen Dif f erentialgleichungssy- 
stem "unterscheidet", wie viele Ableitungsschritte erforder- 
lich sind, um aus dem dif f erential-algebraischen Gleichungs- 
system ein explizites gewdhnliches Dif f erentialgleichungssy- 
stem zu erhalten. 

Ohne Einschrankung der Allgemeingtiltigkeit und bei Existenz 
verschiedener Begrif fsdef initionen des Begriffs "Index", wird 
im weiteren folgende Definition fur den Index eines differen- 
tial-algebraischen Gleichungssystems verwendet: 

Gegeben sei ein dif ferential-algebraisches Gleichungssystem 
des Typs 



20 



F x(t), x(t), t 



= 0 



(2) 



Existiert eine kleinste nattlrliche Zahl i, so dali die Glei- 
chungen 



x(t), x(t), t 



= 0, 



(5) 



25 



lf|x(t), x(t), tj 



dt 



= 0, 



(6) 



d^Bj x(t), x(t), t 



30 



dt 



(i) 



= 0 



(7) 
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algebraisch in ein System expliziter gewohnlicher Differenti 
algleichungen umgeformt werden konnen, so wird i als der In- 
dex des dif ferential-algebraischen Gleichungssystems bezeich 
5 net. Die Funktion F wird dabei als ausreichend oft differen- 
zierbar vorausgesetzt . 

Unter einem stochastischen Dif f erentialgleichungssystem ist 
im weiteren allgemein folgendes Dif ferentialgleichungssystem 
10 zu verstehen: 

Gegeben sei ein Wiener-Hopf-Prozefi jw t ; t e 91q| auf einem 

Wahrscheinlichkeitsraum (Q, A, P) zusammen mit seiner kanoni- 
schen Filtration <{c s ; s e [a, b]} . Seien f erner. h und G:. 

15 [a,b] x 9t -» 5R zwei (( B [ a ,b] x B ) " b) -meflbare Zuf all'svariablen 
und X:f2 -» <R eine (c a - B)-meIibare Funktion. Ein stochasti- 
sches Dif ferentialgleichungssystem ist gegeben durch, das Ito 
Differential- . • - 



25 



20 X s : n -» 9* , 



<d i-» x((d) + f h(u, X u (ffl))dX(uJ'+ j" g(u, X u )dW u , 



(©), (8) 



[a,s] 

{ * i 3 

II. 

It 

Oder symbolisch • " 

dX s = h(s, X s )ds + G(s, X s )dW s , s e 5 [a, b]; :: X a = X . : (9) 

Folgendes Verfahren zur^Rauschs'imulation ist aus [2] bekannt 

30 Die Vorschrift (4 J ;ka T nn filr den Fall einer rein additiven 

Storung entkoppelt, werden in einen. dif f erentiellen und einen 
algebraischen Teil. 0 y 

; * 

' T * 

Fur den Fall einer rein ,_ addi t i<ven~ S t o rung gilt: - ,J 
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B(t,x(t)) ■ B(t), (10) 

d.h. die Intensitatsmatrix hangt nur von der Zeit t ab. 
Damit wird (4) zu 

C ■ x(t) + G ■ x(t) + s(t) + B(t) • v(g>, t) = 0 . (ll) 



10 Unter den Annahmen, dafl konsistente Anfangswerte 

zu einer Anfangszeit to gegeben sind und 



Xdetfco)' x det(to)] 



bei Regularitat des Matrixbuschels {XC + G; X € C}, exi- 
stiert eine eindeutige Losung x^et zu (1) in der Form 



15 C • xdet+ G • x det + s = 0 . (12) 

Ein Matrixbuschel {A.C + G; X e C} ist regular, falls ein Xq 
aus C existiert derart, dafi gilt: 

20 det(X 0 C + G) * 0 . (13) 



Konsistente Anfangswerte 



x det( fc o)' x det(to) konnen dadurch 



gewonnen werden, dafi man einen DOArbeitspunkt des Systems, 

welches durch (12) beschrieben wird, bestimmt, also xdet = 0 
25 setzt. Aus [3] ist ein weiteres Verfahren zur Ermittlung kon- 



sistenter Anfangswerte 



x det(to)/ xdet (to) 



bekannt , 



Nach einer Transformation, die im folgenden beschrieben wird, 
gelangt man zu Vorschrif ten, die zur Vorschrift (11) aquiva- 
30 lent sind, der folgenden Form: 
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y"+ F X • yM + F 2 • yM + d 1 ] + (§(t) • v(<d, t))M = 0 



(14) 



und 



5 F 3 • yM + F 4 • yW + J 2 ) + (§(t) • v(a>, t))^ = 0 



(15) 



mit trans formier ten Anfangsbedingungen 



y^ldetCto): = (t" 1 • x det (t 0 )) [l l 



10 



y det(to> = T • xdet(to) 



(16) 



(17) 



und 



15 yt 2 ldet(to): = (t -1 ■ x det (t 0 )) , 



(18) 



•I 2 ) f . • ^[2] 

y det(to) = T 1 • x d et(to) 



(19) 



und mit 



20 



B(t): = S • B(t) . 



(20) 



25 



Mit Fj_ , i = 1, 2, 3, 4, wird jeweils eine vorgebbare Matrix 
bezeichnet . 

Die Transformation hat zum Ziel, die Vorschrift (11) in ein 
semi-explizites Dif ferentialgleichungssystem des Typs der 

Vorschriften (14) und (15) in der Variable y = (yN, y[ 2 ij mit 
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geeigneten Matrizen F^ und einer Funktion a zu 
ttberfiihren. 



Zur Matrix C aus (11) findet man dazu zwei regulare Matrizen 
5 S und T derart, dafl die Vorschrift 



S • C • T = blockdiag(l, N) 



(21) 



erfUllt ist, wobei mit I eine Einheitsmatrix der Dimension r 
10 und mit N eine Nullmatrix der Dimension (n-r) bezeichnet 
wird. 

Mittels GauJS Elimination mit vollstandiger Pivotstrategie 
werden zwei regulare Matrizen und Qi ermittelt mit 



15 



20 



25 



PL • C • Qi = C lt 



(22) 



wobei eine Matrix eine rechte obere Dreicksmatrix ist, 

deren r erste Diagonalelemente ungleich dem Wert Null sind. 
Die Matrix hat ab der (r+l)-ten Zeile einschliefilich nur 
Eintrage mit dem Wert Null. Die Matrix Qi wird als orthogo- 
nale Spalten-Permutationsmatrix gewahlt. Die Matrix P^ ist 
das Produkt aus einer linken unteren Dreiecksmatrix und einer 
orthogonalen Zeilen-Permutationsmatrix . 

Durch eine Multiplikation mit einer regularen rechten oberen 
Dreiecksmatrix M^ von rechts werden alle Eintrage der Matrix 
Ci oberhalb ihrer Diagonalen eliminiert: 



30 C2: = Ci • Mi = diag 



(n - rj-xaal) 



(23) 



Durch eine Multiplikation mit einer regularen Diagonalmatrix 
M2 von links werden alle nichtverschwindenden Diagonalele- 
mente der Matrix C2 auf den Wert 1 transf ormiert : 
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C3: = M2 • C2 = diag 



f ^ 

X/ ••■ /l / 0/ ••• /0 
^r-mal (n-r)-mal> 



(24) 



Die Matrizen S: = M2 • und T: = Qi • Mi leisten das Gewiinsch- 
5 te. 



Durch Setzen von y: = T 1 • x , E: = C3 = S • C • T, F: = S • G • T 
und o = S • s in Vorschrift (11) ergibt sich 



10 E-y+F-y + a + B- v = 0. (25) 

Um die spezielle Struktur der Matrix C3 auszunutzen, wird y 
in einen ersten Vektor yW, der die ersten r Komponenten 
enthalt, und in einen zweiten Vektor y l ^ J , der die ■ restlichen 
15 (n-r) Eintrage enthalt, unterteilt: 

y = (yWy[2]) T . (26) 

Die Matrix F wird in 4 Untermatrizen Fj_ , i = 1, 2, 3, 4 der 
20 Dimensionen r x r, r x (n-r), (n-r) x r, (n-r) x (n-r) aufge- 
teilt : 




(27) 



25 Fur die Matrix E wird eine entsprechende Aufteilung gewahlt. 

Die Matrix ist eine Einheitsmatrix der Dimension r x r 
und die Matrizen E2, E3 und E4 sind Nullmatrizen. Somit zer- 

fallt (25) in die beiden Vorschriften (14) und (15). 



30 
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1st die Matrix F4 invertierbar, was genau dann der Fall ist, 
wenn das System aus Vorschrift (12) den Index 1 besitzt, so 
kann die Vorschrift (15), 

5 F 3 • yW + F 4 • y[ 2 l + J 2 1 + (§(t) • v(©, t))t 2 J = 0 (15) 

nach yl 2 ] aufgelost werden, was zu folgender Vorschrift 
f Uhrt : 



10 yl^J = -F 4 



F 3 • y M + J2] + (g(t) • v((D, t))t 2 i} . (28) 



Im folgenden werden folgende abkiirzende Bezeichnungen einge- 
f Uhrt : 

15 F: = Fi - F 2 • F4 _1 ■ F 3 ; (29) 
s:= J 1 ! - F 2 • F 4 _1 ■ al 2 l; (30) 
B(t, x): = (l r ,-F 2 • F4- 1 )" 1 • §(t) . (31) 

20 

Mit I r wird eine Einheitsmatrix der Dimension r, also des 
Rangs der Matrix C bezeichnet. 

Durch Einsetzen von Vorschrift (28) in Vorschrift .(14) ergibt 
25 sich: 

M - hi 

y + F • yl l l + s + B(t) • v(<B, t) = 0 . (32) 

Vorschrift (32) kann als ein stochastisches Dif f erentialglei- 
30 chungssystem interpretiert werden der folgenden Form: 

dyf 1 ! = yW - |f • yJ 1 ! + sjdt - B(t)dW t . (33) 
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Mit Y£ J wird eine Zuf allsvariable mit einem Erwartungswert 

y[^| t (to) bezeichnet, die eine endliche Varianz aufweist. 
|v?t;t € 5Rq| ist ein Wiener-Hopf-Prozefi der Dimension der An- 
zahl der Rauschquellen, allgemein der Anzahl der Storquellen. 

Das Verfahren aus [2] geht aus von Vorschrift (33), deren 
eindeutiger Losungsprozefi y£^ als Ito-Dif f erential durch die 
Gleichung 



io = * t#t() • yW - 



J <D t/U • s(u)dX(u) + J <D t/U • B(u)dW u 



■■ v[ fc 0't] . t 0 . „ 

■r 

mit dem Fundamentalsystem von Losungen 



(34) 



15 



25 



*t,t 0 = ex P 



j Fd^(u) 

; [t 0 ,t] . 



= exp{(t 0 - t)F} 



.. ,T35) 



gegeben ist. 



Bei dem Verfahren aus [2] Werden die Erwarttihgswerte Ej- und 
die zweiten Momente der ZufallsvarialDlen naherungs- 

20 weise bestimmt. Der Erwartungswert eines* rtfo-Integrals ist 
gleich dem' Wert 0. Somit erhalt m'an aus.,,Vorschrif t (34)/di- 
rekt 



E t := e(y[ 11 ) = *t,t 0 • yStfro) " \ *t,u • i(u)dX(u). 



(35) 



Fur alle t lost Et' das gewohnl'iche Dif f erentialgleichungssy- 
stem 



t r 
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xt+ F • xt + s = 0; 



11 



x t 0 



= y [ dit(to) • 



(36) 



Die zweiten Momente 



5 Pt = E 



(37) 



10 



der Zufallsvariablen des Losungsprozesses gentigen fur 
alle t dem Dif ferentialgleichungssystem 



3j = -F • P t - P t • F T - s • Ej - E t • s T + B(t) • B(t) T , 



:38) 



15 



wobei eine Anf angsbedingung durch 



p tn = E 



i 1 • 0$ 



(39) 



gegeben ist. Bei (38) handelt es sich urn ein lineares gewohn- 
liches Dif ferentialgleichungssystem. 



Parallel zur transienten Simulation der Schaltung werden bei 
20 dem Verfahren aus [2] durch numerische Integration mit linea- 
ren impliziten Mehrschrittverf ahren naherungsweise die Erwar- 
tungswerte und die zweiten Momente P^- ermittelt. 



Ein Nachteil dieses Verfahrens ist darin zu sehen, dafi zur 
25 Bestimmung der zweiten Momente Pt fiir jeden Zeitschritt li- 
neare Gleichungssysteme von quadratischer Ordnung in der An 
zahl m der Rauschquellen gelost werden. 



Das Verfahren basiert auf einer manuellen Indexreduktion des 
30 dif f erential-algebraischen Gleichungssystems auf ein explizi- 
tes stochastisches Dif ferentialgleichungssystem, die nicht 
automatisierbar ist. AuJlerdem wird nur ein Teil der Rauschef- 
fekte betrachtet. Ferner liefert das Verfahren keine pfadwei- 
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se Information, sondern nur die Momente der Knotenpotentiale, 
die bei der Indexreduktion erhalten bleiben. Fur die Knoten- 
potentiale, die durch die Indexreduktion unterdruckt werden, 
liefert dieses Verfahren keine Information. 

5 

Ein weiterer Nachteil des Verfahrens aus [2] ist darin zu se- 
hen, dafi bei diesem Verfahren die Indexreduktion aulierst 
ineffizient sowie auch nicht automatisch durchfuhrbar ist, 
insbesondere da die algebraischen Variablen in dem differen- 
10 tial-algebraischen Gleichungs system nicht beriicksichtigt wer- 
den. Die Indexreduktion mull bei dem Verfahren aus [2] analy- 
tisch manuell durchgefiihrt werden, da die numerischen Verfah- 
ren nicht stabil sind. 

15 Ferner ist es aus [4] bekannt, eine Rauschsimulation einer 
integrierten Schaltung im Frequenzbereich durchzufuhren. Da- 
durch lalit sich eine Schaltung jedoch nur im Kleinsignalbe- 
reich analysieren. Die Voraussetzung eines festen Ar- 
beitspunktes ist jedoch haufig nicht gegeben. Beispielsweise 

20 verhindert das Schwingverhalten eines Oszillators in einer 
Schaltung einen einheitlichen Arbeitpunkt. 

In [5] ist eine Schaltungsbeschreibungssprache SPICE be- 
schrieben, mit der eine elektrische Schaltung in einer ftir 
25 einen Rechner verarbeitbarer Form beschreibbar ist. 

Aus [6] ist ein Verfahren zur numerischen Behandlung eines 
stochastischen Dif f erentialgleichungssystems bekannt, die 
pfadweise Simulation diskreter Approximationen an den L6- 
30 sungsprozefi, das als Runge-Kutta-Schema bezeichnet wird. 

Aus [7] ist ein Simulator bekannt, der 
o eine Initialisierungseinheit, 
o eine Inkrementiereinheit, 
35 o eine Einheit zur Aktualisierung eines Schatzwerts, und 
o eine Ausgabeeinheit 
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umfaiit. Von der Initialisierungseinheit wird ein Initialisie- 
rungswert eines Zustands und einer Funktion, mit denen ein 
System, welches einer zufalligen Storung unterliegt, be- 
schrieben wird, der Inkrementiereinheit zugefUhrt. Die Inkre- 
mentiereinheit verwendet zwei Zuf allszahlenf olgen urn ein In- 
krement des Schatzwerts zu bilden, ohne die Funktion selbst 
zu dif f erenzieren. In der Einheit zur Aktualisierung eines 
Schatzwerts wird der Schatzwert urn das Inkrement erhoht. 

Weiterhin sind aus [8] ein Verfahren sowie eine Vorrichtung 
zur Verbesserung der Genauigkeit eines in einem geschlossenen 
Regelkreis geregelten Systems bekannt, wobei mindestens zwei 
stochastische Storsignale beerticksichtigt werden. 

Somit liegt der Erfindung das Problem zugrunde, ein Verfahren 
anzugeben, mit dem die im vorigen beschriebenen Nachteile 
vermieden werden. 

Das Problem wird durch das Verfahren gemali Patentanspruch 1 
und durch die Vorrichtung gemali Patentanspruch 12 gelost. 

Bei dem Verfahren gemafi Anspruch 1 wird ein technisches Sy- 
stem, welches einer StSrung unterliegt, mittels eines impli- 
ziten stochastischen Dif f erentialgleichungssystems (SDE) be- 
schrieben. Eine naherungsweise Losung des Systems wird ermit- 
telt, indem ein diskreter NaherungsprozeJi realisiert wird. 
Der diskrete Naherungsprozeii wird gemafi folgender Vorschrift 
realisiert: 




+B(x n ) • AW n + yB(x n ) • AWn_i 



wobei mit 



— C eine erste Matrix, 
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— <xi, ct2/ Y vorgegebene Parameter aus dem Intervall [0, 1], 

T 

— h: = — , eine Schrittweite in emem Ausgangszeitintervall 

N 

[0, T], wobei T ein vorgegebener Wert ist, welches in N 
Teilintervalle unterteilt ist, 
5 — G eine zweite Matrix, 

— ^x n +i eine Realisierung des Naherungsprozesses an einer 

Stutzstelle x n+ i, 

— X Tn eine Realisierung des Naherungsprozesses an einer 

Sttltzstelle x n , 

10 — *x n _i eine Realisierun< ? des Naherungsprozesses an einer 
Stiitzstelle x n -l' 

— s(t n + i) ein erster Wert an der Stutzstelle x n +l' 

— s(x n ) ein erster Wert an der Sttltzstelle x n , 

— s(x n _i) ein erster Wert an der Stutzstelle x n -i, 

1,5 — ein Dif f erenzwert = W Tn - W x ^ zwischen einem zwei- 

ten Wert W tn an der Stutzstelle x n und dem zweiten Wert 

W T n an der Stutzstelle x n _! 
c n — 1 

— B(x n ) ein zweiter Wert an der Stutzstelle x n/ 

bezeichnet wird, und 
20 - bei dem durch iterative Losung des Naherungsprozesses die 
Storung X Xn+1 ermittelt wird. 

Die Vorrichtung gemafi Patentanspruch 12 weist eine Prozes- 
soreinheit auf, die derart eingerichtet ist, dafi ein techni- 

25 sches System, welches einer Storung unterliegt, mittels eines 
impliziten stochastischen Dif f erentialgleichungssystems be- 
schrieben werden. Eine naherungsweise Losung des Systems wird 
ermittelt, indem ein diskreter Naherungsprozell realisiert 
wird. Der diskrete Naherungsprozell wird gemafl folgender Vor- 

30 schrift realisiert: 
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— (c + ha 2 c) • X Tn + 1 = {- (l - y)c + h(l + yet! - a 2 )c} ■ X Xn + 

+{- yC + h(y - aiyjc} . X^^ + 
+ha 2 s(x n + 1 ) + h(l + yai - a 2 )s(x n ) + 
+h(y - aiY)s(x n „!) + 
+B(x n ) • AWn + YB(x n ) • A%_i 

wobei mit 

— C eine erste Matrix, 

— ai, a 2/ y vorgegebene Parameter aus dem Intervall [0, 1], 

T 

5 — h: = — , eine Schrittweite in emem Ausgangszeitintervall 

[0, T], wobei T ein vorgegebener Wert ist, welches in N 
Teilintervalle unterteilt ist, 

— G eine zweite. Matrix, 

— ^x n +i eine Realisierung des Naherungsprozesses an einer 

10 Stiitzstelle x n+ i, 

— X T eine Realisierung des Naherungsprozesses an einer 

c n 

Stiitzstelle" x n , 

— X T «, eine Realisierung des Naherungsprozesses an, "einer 

Stiitzstelle x n _i, 
15 — s(x n + i) ein erster Wert an der Stiitzstelle x n + l' 

— s(x n ) ein erster Wert an der Stiitzstelle x n , ... - 

— s(x n _i) ein erster Wert" an der Stiitzstelle x n _i, 

— ein Dif ferenzwert A%1].: = W Xn - w x n _i zwischen einem zwei 

ten Wert W T an der Stiitzstelle x n und-dem zweiten Wert 
n 

20 W x . an der Stiitzstelle x n _i . . *v 

^n-l 11 < /• " 

— B(x n ) ein zweiter Wert an der Stiitzstelle x n , ' 
bezeichnet wird, Durch iterative*' Losung des Naherungsprbzes 
ses wird die Storung Xx n +i ermittelt. 

25 Die Erfindung verwendet direkt die implizite Struktur des 
technischen Systems, represent ier't durch ein implizites Dif 
ferenti^JLgleichuHgssystem. DUtch'die Erfindung wird die Er- 
mittlung der Storung erheblich beschleunigt, da die Dunnbe- 
setztheitc der Matrizen C und" J G ausgenutzt werden kann. ,pie - 



I u 

f 0 " 
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numerisch instabile und aufwendige Transformation des Diffe- 
rentialgleichungssystems in die entkoppelte Form entfallt. 

Durch die Erfindung ist es erstmals moglich, eine Storung 
5 auch bei einer singularen Matrix C zu ermitteln. 

Vorteilhafte Weiterbildungen der Erfindung ergeben sich aus 
den abhangigen Anspriichen. 

10 Die Erfindung kann in verschiedensten Anwendungsgebieten ein 
gesetzt werden, in denen ein technisches System gestort ist 
und durch ein System dif f erential-algebraischer Gleichungen 
beschrieben werden kann. 

15 Beispielsweise konnen Storungen (Rauschen) in einer elektri- 
schen Schaltung ermittelt werden. Die Erfindung eignet sich 
ferner zur Anwendung in einem mechanischen Mehrkdrpersystem 
oder in einem allgemeinen physikalischen System, einem chemi 
schen System oder auch einem physikalisch-chemischen System, 

20 deren jeweilige Modellierung auf ein System dif f erential- 
algebraischer Gleichungen ftihrt. 

Ohne Einschrankung der Allgemeingiiltigkeit wird im weiteren 
die Erfindung anhand eines AusfUhrungsbeispiels einer Rausch 
25 simulation einer elektrischen Schaltung beschrieben, das in 
den Figuren dargestellt ist. 



Es zeigen 

Figur 1 ein Ablaufdiagramm, in dem die einzelnen 
30 Verfahrensschritte dargestellt sind; 

Figur 2 eine Vorrichtung, mit der das Verfahren durchge- 
fuhrt wird; 

Figur 3 eine Skizze einer elektrischen Schaltung eines 
Differentiators; 
35 Figur 4 eine Darstellung eines Simulationsergebnisses . 
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In Figur 2 ist in Form einer Blockskizze ein Rechner R darge- 
stellt, mit dem das im weiteren beschriebene Verfahren durch- 
gefUhrt wird. Der Rechner R weist eine Prozessoreinheit P 
auf, die derart eingerichtet ist, dali die im weiteren be- 
5 schriebenen Verf ahrensschritte durchgef Uhrt werden konnen. 
Die Prozessoreinheit P ist uber einen Bus B mit einem Spei^ 
cher SP verbunden. Der Rechner R ist uber eine erste Schnitt- 
stelle I/O 1 mit einer Tastatur TA und uber eine zweite 
Schnittstelle I/O 2 mit einer Maus MA verbunden, mit denen 
10 jeweils ein Benutzer des Rechners R Eingaben vornehmen kann. 
Ferner ist der Rechner R mit einem Bildschirm BS verbunden, 
auf dem dem Benutzer Ergebnisse des Verfahrens dargestellt 
werden. 

15 In dem Speicher SP ist die zu analysierende Schaltung S 

(vgl. Figur 3) in Form einer Schaltungsbeschreibungssprache 
gespeichert. Als Schaltungsbeschreibungssprache wird das 
Netzlistenformat von SPICE eingesetzt. 

20 Ein Ersatzschaltbild der Schaltung S fttr dieses Ausfuhrungs- 
beispiel ist in Figur 3 dargestellt. Die Schaltung S stellt 
in ihrer Funktionalitat einen Differentiator dar . Eine unab- 
hangige Spannungsquelle V(t) ist zwischen einem ersten Knoten 
und einem Massepotential Nq angeordnet. Zwischen dem er- 

25 sten Knoten und einem zweiten Knoten N2 ist eine Kapazi- 
tat Cq vorgesehen. Der zweite Knoten N2 ist uber einen Wi- 
derstand R mit einem dritten Knoten N3 verbunden. Ferner ist 
zur Modellierung eines thermischen Rauschens in dem Wider- 
stand R eine zu dem Widerstand R parallel geschaltete Rausch- 

30 stromquelle RS mit der Stromstarke Al R vorgesehen. Ein er- 
ster Eingang El eines Operationsverstarkers OV ist mit dem 
zweiten Knoten N2 verbunden, ein zweiter Eingang E2 des Ope- 
rationsverstarkers OV mit dem Massepotential Nq . Ein Ausgang 
A des Operationsverstarkers OV ist mit dem dritten Knoten N3 

35 verbunden. 
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Die Schaltung S wird in einem ersten Schritt 101 in Form des 
Netzlistenformats von SPICE in dem Speicher SP gespeichert 
(vgl. Figur 1) . 

5 In einem zweiten Schritt 102 wird filr die Schaltung S-" eine 
modifizierte Knotenanalyse durchgefuhrt . Das Ergebnis der mo- 
difizierten Knotenanalyse ist das der Schaltung S zugehorige 
Gleichungssystem. 



10 In Matrixschreibweise ergibt sich fur die Schaltung S das li- 
near-implizite Dif f erentialgleichungssystem des Indexes 1 mit 
rein additiver Storung 



15 



20 
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Es werden mit 

- ui, U£/ U3 die Spannungen an den Knoten N]_, N2, N3 , 

- Iq der durch den ersten Knoten und durch die 

Spannungsquelle V(t) flieliende Strom, 

- I^pi der durch den dritten Knoten N3 und den 

Operationsverstarker OV flieliende Strom, 

- R der Wert des Widerstand R, 

- C 0 der Wert der Kapazitat C 0 , 
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- A ein Verstarkungsf aktor des Operationsverstarkers OV, 
bezeichnet . 

In dem Modell des thermischen Rauschens ist die Stromstarke 
5 AIr des Stroms,;der durch die Rauschstromquelle RS flielit, 

durch die Gleichung 



20 



/ 4k T 

AI R = J— Af ■ v(cd, t) (40) 

10 angenommen . 

Hierbei beschreibt v(<d, t) als weilies Rauschen den treibenden 
verallgemeinerten stochastischen ProzeJJ, also die Storung des 
technischen Systems. Die Grofien k, T und Af sind als Kon- 
15 stanten vorgegeben. 

Allgemein ergibt sich also ein nichtentkoppeltes linear- 
implizites Dif ferentialgleichungssystem mit rein additiver 
stochastischen Storung (Rauschen) der folgenden Form: 



C ■ x(t) + G ■ x(t) +' s(t) + B(t) ■ v(©, t) = 0 . (41) 



Zur Herleitung des Verfahrens, nicht aber fur das Verfahren 
selbst, wird angenommen , daJ3 die Kapazitatsmatrix C inver- 
25 tierbar ist. 

Durch Multiplikation der Gleichung (41) mit der inversen Ka- 
pazitatsmatrix C" 1 erhalt man 

30 x(t) = -CT 1 • {G • x(t) + s(t)} - CT 1 ■ B(t, x(t)) ■ v(<D, t) . (42) 

Die Vorschrift (42) wird als stochastisches Dif ferentialglei- 
chungssystem der Form 

35 dX t = X tQ - CT 1 • {G • X t + s}dt - CT 1 • B(t)dW t (43) 
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interpretiert, wobei Xt Q eine Zuf allsvariable mit Erwar- 
tungswert Xd e t(to) bezeichnet. 

5 Auf das Dif ferentialgleichungssystem (37) wird folgendes Ver- 
fahren, das in [6] beschrieben ist, angewendet (Das Verfahren 
wird als implizites starkes Zweischrittverf ahren der Ordnung 
1 ftir ein stochastisches Dif ferentialgleichungssystem be- 
zeichnet) : 

10 

Ein numerisches Verfahren liefert fur jede naturliche Zahl N 
Realisierungen eines diskreten Naherungsprozesses 
{x s ; s = xo/..-,x n }, indem es fur gegebene cd aus Q an den je- 
weiligen Stiitzstellen x-^ Approximationen x(co, xjj fur die Wer- 
15 te x(o, XjJ berechnet. 

Um einen stetigen Naherungsprozefi {x s ;s <= [0, T]} zu erhalten, 

interpoliert man die erhaltenen Werte z.B. stuckweise linear. 
Die sich ergebenden Pfade konnen anschliefiend mit statisti- 
20 schen Methoden analysiert werden. 

In [6] ist folgende Familie von impliziten starken Zwei- 
schrittverf ahren der Ordnung 1 zur Bestimmung der Komponente 
k der Realisierungen |x Ti ; i = 0, ...,nJ angegeben: 

25 

*n + l = (l - Vk)*£ + Yk^ n ^ 1 

^kf^n+l^ X Tn+1 ) + (l + Yk<*l,k * a 2/k )f k 

+ Tk(l " <*1, k)f ^n-1^ Xtn-i) 
+ v £ + YkV^-i 

Hierbei steht eine Funktion ohne Argumente fur eine Auswer- 
tung am Punkt ( T n' *x n ) • D:i - e reellen Parameter ai, k, ot2,k und 

30 7k sind aus dem Intervall [0, 1] gewahlt. Ferner ist die Gro- 
k 

]Je V n gegeben durch 
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3 = 1 31.32 

wobei 

AW n := W Xn + 1 - W Xn (46) 

die Zuwachse des Wiener-Prozesses beschreibt und I? eine ab- 
kiirzende Schreibweise fur den Operator 

d . d 

!?:= 2]g k ' j —xt, j = 1, ...,m (47) 



k=l 5x1 ( 

' - . ! J 1 « 

ist. 

15 Die Grolie I/. . \ beschreibt das mehrfache Ito- 

gi/J2j' x n' x n + l 

Integral - . - 

i i M , . ; / v * 

' i 

1 • ■ . 

*n + l s 2 . . ' 
If. . \ := f f dwjdwj . (48) 

j * \ ■ ' 

t ■ ..'to. 

20 h ist gegeben durch die Vorschrift 

- » 

T J v : 1 

h:= - . ,. ' « .; " T49) 

N , . ' 

ft • ' / 1 - * ■ • 

Wahlt man fur alle k = 1; ... , n die Parameter ai,^, ot2,k und 
25 Yk, so vereinfacht siqh Vorschxift (44) fiir den Fall rein ad- 
ditiver Storung in Vektorschreibweise zu 

+h{a 2 f(x n +l/Xx n+1 ) + (1 + yai - a 2 )f + : y(l - a^in-i, X^')} , 

+g • AW n + yg - AW n _i___ " .. ^ 

- ■** * — - — - — _ • * * r, v 

■ ji ■ - Iff ■ ■ _ | ft 



WO 99/48030 



PCT/DE99/00239 



22 



(50) 



-h 



da die Funktion g nur von der Zeit abhangt. Vorschrift (50) 
wird auf das Dif f erentialgleichungssystem (43) angewendet. 

Es ergibt sich somit: 

x Xn+1 = (i - y^x n * yxx^! 

<x 2 C~ 1 (g • X Tn + 1 + s(x n+ i)) + 
+(l + yai - a 2 )c _1 (G • X Xn + s(x n )) . (51) 
+(y - axyJc-^G • X Xn _ 1 + s(x n _i)) 
-C" 1 B(x n )AW n - yC-^TnjAW^-i 

10 Im nachsten Schritt wird (40) mit -C multipliziert . 
Es ergibt sich somit: 

- c K + i -( l ~y)K "An-lK 

a 2{ G 1 Xx n+i + s(x n + 1 )} 

+ (l + yai - cx 2 ){g • X Xn + s(x n )} , (52) 
+ (y - a lY ){G • X Xn _ 1 + s(x n _i)[ 
+B(x n ) • AWn + yB(x n ) • AW n _! 

15 

Durch Zusammenfassen der Terme in X» auf der rechten Sei- 

k n+i 

te von (41) ergibt sich: 
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- (c + ha 2 c) • X Xn + 1 = {- (l - y)c + h(l + yai - a 2 )c} ■ X Tn + . 

+{- yC + h(y - aiy)c} • X Xn ^ 1 + 
+ha 2 s(x n + i) + h(l + yai - a 2 )s(x n ) + 
+h(y - a 1 y)s(x n . 1 ) + J 
+B(x n ) • A% + yB(x n ) ■ AW n _ x 

(53) 

Ftir den Fall des Indexes 0 ist durch Vorschrift (53) eine 
5 Vorschrift angegeben, mit der die Storung des Systems ermit- 
telt werden kann. Vorschrift (53) kann auch im Fall des Inde- 
xes 1 angewendet werden. 

Durch iterative Losung des Naherungsprozesses wird in einem 
10 weiteren Schritt 103 die Storung Xt n+1 ermittelt. Die Be- 

stimmung der Storung Xt n + 1 wird iterativ in einer Schleife 

flir n = 0,1/ ... , N - 1 durchgeftihrt, wodurch der Verlauf der 
Storung X t j zu den jeweiligen Zeitpunkten x n ermittelt 

wird. 

15 

Abhangig von der ermittelten Storung wird eine Schaltung S 
modifiziert, so dafi die vorgebbaren Bedingungen z.B. hin- 
sichtlich der Storanf alligkeit der Schaltung S erfullt sind. 

20 Die modif izierte Schaltung wird in einem letzten Schritt in 
Silizium gebrannt. 

Figur 4 zeigt als Ergebnis des von dem Rechner R durchgefiihr- 
ten Verfahrens, das im weiteren als Programm angegeben ist. 
25 Es ist dargestellt ein numerisch simulierter Losungspfad RLP 
des Spannungsverlaufs im dritten Knoten N3 unter Berucksich- 

tigung des Rauschens. Ferner ist zum Vergleich der ideale Lo- 
sungspfad ILP, d.h. der Losungspfad ohne Berucksichtigung des 
Rauschens angegeben. AuJierdem ist der der Verlauf V der Ein- 
30 gangsgrdfle V(t) dargestellt. 
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Fur das Verfahren wurden folgende Parameter verwendet: 



A = 2,5 • 10 



V(t) = 



2 • io c 



-2 • 10 



C 0 = 1 
R = 1 - 10 
Af = 1 



10 
4 



-12 



0 far 0 < t <, 5 • 10' 



-9 



t - 1 
,-2 



10 



-2 



fur 5 • 10 9 < t < 10 • 10~ 9 



,-9 



,-9 



10 * ftir 10 • 10 ' < t < 15 • 10" 
t + 4 • 10" 2 fur 15 • 10" 9 < t < 20 - 10~ 9 
0 



filr 20 - 10 9 < t < 25 • 10" 9 



[s] 
[V] 



[F] 
[Hz] 



Im weiteren sind einige Varianten bzw. Verallgemeinerungen 
des oben beschriebenen Ausftlhrungsbeispiels dargestellt: 



Die Kapazitatsmatrix kann ohne weiteres auch singular sein. 
Somit ist es erstmals mSglich, auch eine singulare Kapazi- 
tatsmatrix zu berucksichtigen. 

10 

Ferner kennen auch weitere in [6] beschriebene Verfahren zur 
pfadweisen Approximation eines stochastischen Dif f erential- 
gleichungssystems eingesetzt werden, urn z.B. eine hohere Kon- 
vergenzordnung zu erreichen. 

15 

Die Erfindung ist keineswegs auf die Ermittlung thermischen 
Rauschens beschrankt. Im Gegenteil kann jedes Modell einer 
Storung im Rahmen der Erfindung berticksichtigt werden, wel- 
ches durch weifies Rauschen beschrieben oder approximiert wer- 
20 den kann, z.B. das Schrotrauschen in Stromquellen. Bei Be- 
rticksichtigung des Schrotrauschens wird die Stromquelle mo- 

delliert durch eine parallel geschaltete, zuf allsabhangige 
Stromquelle, deren Stromstarke AIq der Vorschrift 



25 AIq = ^eldetAf • v(oo, t) 
gentigt . 
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Im weiteren ist eine Realisierung des Verfahrens zur Rausch- 
simulation in der Programmiersprache FORTRAN7 7 angegeben. 



PROGRAM integr 
5 C 

C Verfahren zur numerischen Loesung linear-impliziter 

C Dif ferentialgleichungssysteme mit rein additiver Storung. 

C 

IMPLICIT NONE 
10 INTRINSIC dble, real 

EXTERNAL setall, gennor, dgesv, giveB, gives 
EXTERNAL a_smsm, m_mv, v_to_v, m_sv, a_vvvv 
REAL gennor 

C 

15 C Konstanten- und Variablendeklarationen und Initialisierung 
C 

INTEGER n, m, anz, i, j, info, initl, init2 

C 

C n ist die Dimension des Systems, m die Anzahl der 
20 C Rauschquellen 

C anz ist die Anzahl der betrachteten Teilintervalle, 

C i und j sind Hilf svariablen fur Schleifen, 

C info fragt den Return-Code 

C initl und init2 sind die Samen fuer den 
25 C Pseudozufallszahlengenerator 

C 

PARAMETER (n = 5, m = 1, anz = 1000) 
PARAMETER (initl = 23556, init2 = 4285979) 

C 

30 DOUBLE PRECISION Irand, rrand, h, tauakt, alpha 

C 

C [Irand, rrand] bezeichnet das betrachtete Intervall 
C h die Schrittweite, tauakt den aktuellen Zeitpunkt, 
C alpha den Parameterwert alpha2 aus (53), wobei alphal = 0.0 
35 C und gamma =0.0 gesetzt wurde. 
C 

PARAMETER (alpha =0.9) 
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PARAMETER (lrand ■ O.OdO, rrand = 2.5d-8) 

C 

DOUBLE PRECISION xakt (n) , d(n), salt (n) , sakt(n) 
DOUBLE PRECISION C(n,n) , G(n,n), B(n,m), A(n,n) 
5 DOUBLE PRECISION hl_n (n) , h2_n(n), h3_n(n), h4_n(n) 

DOUBLE PRECISION hl_m(m) , hl_nn(n, n) 

C 

C xakt beschreibt den Wert des Losungsprozesses am Zeitpunkt 
C tauakt, 

10 C C, B und G sind die Matrizen des Problems der Vorschrift 
C (3) 

C A und d dienen zum Gleichungsaufbau 
C salt und sakt beschreiben s(taukt) bzw. s(takt - 1) 
C h*_* bezeichnen Hilfsvektoren und -matrizen 
15 C 

C Konsistenten Anfangsvektor setzen 
xakt(l) = O.OdO 
xakt (2) = O.OdO 
xakt (3) = O.OdO 
20 xakt (4) = O.OdO 

xakt (5) = O.OdO 

C 

C Definition der Matrix C (mit Kapazitat = 1* 10 A {-12} 
C Farad) 



25 


C(l, 1) 




1.0d-12 




C(2, 1) 




-1.0d-12 




0(3,1) 




O.OdO 




0(4,1) 




O.OdO 




C(5,l) 




O.OdO 


30 


C(l,2) 




-1.0d-12 




0(2,2) 




1.0d-12 




0(3,2) 




O.OdO 




0(4,2) 




O.OdO 




C(5,2) 




O.OdO 


35 


0(1,3) 




O.OdO 




C(2,3) 




O.OdO 




0(3,3) 




O.OdO 
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C (4, 3) 




0 . Odu 




C (5, 3) 




0 . OdO 




C (1, 4) 




A A -3 A 

0 . OdO 




C (2, 4) 




0 . OdO 


5 


C (3, 4) 




A A -J A 

0 . OdO 




C (4, 4) 




A A J A 

0 • OdO 




C (5, 4) 




A A -J A 

0 . OdO 




C (1/ 5) 




u . uau 




C (2, 5) 




0 . OdQ 


i n 
IU 


C (O, 3) 




O.OdO 




U 1*1/ DJ 




0 . OdO 




^ ^ O/ o j 




O.OdO 












l jjeriniLion 


der Matrix 


ID 


G (1/ 1)" 




[' O.OdO 




G C2; 1) 




0 . OdO 




G (3, 1) 




0 . OdO 




G (4, 1) 




O.OdO 




G (5, 1) 




1 OHO 


20 


G ( 1> 2 ) 

if 




O.OdO 




G (.2, 2) 




1.0d-4 




G (3, 2) 




-1.0d-4 




f** t A A \ 

0(4,2) 




3.0d2 




G ( O/ 2 ) 




O.OdO 


2S 


G ( 1/ J) 




O.OdO 




G (2, 3) 




-1.0d-4 




G (3/ 3i) 




1 . Od-4 




G ( 4 , J ) 




l.OdO 




G ( 0/ J ) 




O.OdO 




G ( If 4 ) 




l.OdO " 




G (z, 4 ) 




O.OdO 

f 1 




G ( J / 4 ) 




O.OdO 




G (4 ,.4) 




0..0d6" 




G(5,4) 




0 . OdO . 
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Gtl/5) 




'o.OdQ 




G{2,5) 




O.OdO 




G(3',5) 




-l.OdO. 
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0(4,5) = O.OdO 
0(5,5) » O.OdO 



C 



5 



C 



Pseudozufallszahlengenerator initialisieren 
CALL setalMinitl, init2) 



C 



C 



Schrittweite berechnen 



C 



h = (rrand - lrand) / anz 



10 



C 



C 



Zeitpunkt und x-Wert initialisieren 



C 



tauakt = lrand 



C 



15 C Ausgabedatei offnen 

OPEN (9/ FILE = f opmitlOOO' , FORM = 1 FORMATTED f ) 
C Anfangszeitpunkt und -wert in die Ausgabedatei schreiben 

WRITE (9, 42) tauakt, (xakt(i), i = 1, n) 
C "alten" s-Wert merken 
20 CALL gives (salt, n, tauakt) 

C 

C Aussere Schleife: Entspricht der Behandlung eines 
C Teilintervalls 

DO 10, i = 1, anz 
25 C Aufbauen des Systems A*X=b 
C 

C Aufbauen von A 

C A = -C -alpha * h * G 

CALL a_smsm(C, G, -l.OdO, -alpha * h, A, n, n) 
30 C B aufbauen 

CALL giveB(B, n, m, tauakt) 
C Zeitpunkt erhohen 

tauakt = tauakt + h 
C sakt aufbauen — salt liegt vor Schleif eneintritt vor 
35 CALL gives (sakt, n, tauakt) 

C "alten" s-Wert merken 



CALL v to v(sakt, salt, n) 
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C d aufbauen 
C 

C hl_nn = -C + (1- alpha) * h * G 

CALL a_smsm(C, G, -l.OdO, (1-alpha)* h, hl_nn, n, n) 
5 C hl_n = hll_nn * xakt 

CALL m_mv(hl_nn, xakt, hl_n, n, n) 
C h2_n = alpha * h * sakt 

CALL m_sv(sakt, alpha*h, h2_n, n.) 
C h3_n = (1 - alpha) * h * salt 
10 CALL m_sv(salt, (1- alpha) * h, h3_n, n) 

C hl_m = DeltaW 

DO 20, j = 1, m 

hl_m(j) = DBLE (gennor (0.0, SQRT (REAL (h) ) ) ) 
20 CONTINUE 
15 C h4_n = B * hl_m 

CALL m_mv(B, hl_m, h4_n, n, m) 
C d = hl_n + h2_n + h3_n + h4_n 

CALL a_vvvv(hl_n, h2_n, h3_n, h4_n, d, n) 
C Aufrufen des Gleichungsl"osers (hl_n ist nur Dummy) 
20 CALL dgesv(n, 1, A, n, hl_n, d, n, info) 

C Neuen Wert von xakt setzen 

CALL v_to_v(d, xakt, n) 
C aktuellen Zeitpunkt und x-Wert in die Ausgabedatei schrei 
C ben 

25 WRITE (9, 42) tauakt, (xakt(j), j = 1, 3) 

42 FORMAT(E16.6,E16.6,E16.6,E16.6,E16.6, E16.6) 

10 CONTINUE 

C 

C Ausgabedatei schliessen 
30 CLOSE (9, STATUS = 'keep') 

STOP 
END 

C 



35 



SUBROUTINE giveB(outB, n, m, tauakt) 
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10 



15 



20 



25 



30 



35 



C 
C 

c 



c 
c 



Liefert in outB den Wert der (n x m) -Matrix B 
aus Vorschrift (3) zum Zeitpunkt tauakt 

IMPLICIT NONE 
INTRINSIC dsqrt 

Dummy Argumente 
INTEGER n,m 

DOUBLE PRECISION outB(n, m) , tauakt 

Lokale Variable 

DOUBLE PRECISION k, T, deltaf 
k bezeichnet die Boltzmann-Konstante, T die absolute 
Temperatur, deltaf die Rauschbandbreite 

PARAMETER (k = 1.308d-23, T = 300, deltaf = l.OdO) 
Zeitunabhangige Definition der Matrix B aus Beispiel 3.4 
(mit Rl = 3000, R2 = 4000, R3 = 5000 Ohm 
T = 300 K, k = 1.3807 * 10 A {-23} jK "{-!}, f = 10) 



outB(l,l 

outB(2,l 

outB(3, 1 

outB(4,l 

outB(5, 1 

RETURN 

END 



= O.OdO 

= dsqrt (4*k*T*deltaf/10000.0dO) 
= -dsqrt (4*k*T*deltaf/10000.0dO) 
= O.OdO 
= O.OdO 



c 

c= 

c 

c 
c 
c 
c 
c 



SUBROUTINE gives (outs, n, tauakt) 

Liefert in outs den Wert des (n)-Vektors s 
aus Vorschrift (3) zum Zeitpunkt tauakt 



IMPLICIT NONE 
INTRINSIC sin 
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C 

C Dummy Argumente 
INTEGER n 

DOUBLE PRECISION outs (n) , tauakt 

5 C 

C Lokale Variable 
INTEGER i 

C i beschreibt die aktuelle Zeilenposition von outs 
DO 10, i = 1, n - 1 
10 outs(i) = O.OdO 

10 CONTINUE 

IF (O.OdO .le. tauakt .and. tauakt .le. 5.0d-9) THEN 
outs(n) = O.OdO 
ELSE IF (5.0d-9 .It. tauakt .and. tauakt .le. 10.0d-9) 
15 THEN 

outs(n) = -(2.0d6 * tauakt - 1.0d-2) 
ELSE IF (10.0d-9 .It. tauakt .and. tauakt .le. 15.0d-9) 
THEN 

outs(n) = -(10.0d-3) 
20 ELSE IF (15.0d-9 .It. tauakt .and. tauakt .le. 20.0d-9.) 

THEN 

outs(n) = -(-2.0d6 * tauakt + 4.0d-2) 
ELSE IF (20.0d-9 .It. tauakt .and. tauakt .le. 25.0d-9) 
THEN 

25 outs(n) = O.OdO 

END IF 
RETURN 
END 
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Patentanspruche 

1. Verfahren zur Ermittlung einer Storung eines technischen 
Systems, welches einer Storung unterliegt, 
5 - bei dem das System mittels eines impliziten stochastischen 
Dif ferentialgleichungssystems beschrieben wird, 

— bei dem eine naherungsweise Losung des Systems ermittelt 
wird, indem ein diskreter NaherungsprozeiS realisiert wird, 

— bei dem der diskrete Naherungsprozefi gemaJl folgender Vor- 
10 schrift realisiert wird: 

— (c + ha 2 c) • X Tn + 1 = {- (l - y)c + h(l + yai - a 2 )c} ■ X Xn + 

+{- yC + h(y - aiy)c} • + 

+ha 2 s(x n + 1 ) + h(l + yai - a 2 )s(x n ) + 

+h(y - aiy)s( x n-l) + 

+B(x n ) ; AW n + yB(x n ) • AW n «x 

wobei mit 

— C eine erste Matrix, 

15 — <xi, a 2 , y vorgegebene Parameter aus dem Intervall [0, 1], 
T 

— h: = — , eine Schrittweite in einem Ausgangszeitintervall 

N 

[0, T] , wobei T ein vorgegebener Wert ist, welches in N 
Teilintervalle unterteilt ist, 

— G eine zweite Matrix, 

20 — **n + l e * ne Realisierung de s Naherungsprozesses an einer 
Stutzstelle x n + i, 

— X Tn eine Realisierung des Naherungsprozesses an einer 

Stutzstelle x n , 

— ^x n _i eine Realisierung des Naherungsprozesses an einer 
25 Stutzstelle x n _i, . 

— s(x n+ i) ein erster Wert an der Stutzstelle x n + i, 

— s(x n ) ein erster Wert an der Stutzstelle x n , 

— s(x n _i) ein erster Wert an der Stutzstelle x n -i, 
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— ein Dif ferenzwert AWn_i: = W Xn - W tn _ 1 zwischen einem zwei- 
ten Wert W Xr an cler StUtzstelle T n und dem zweiten Wert 
W Tn-1 an der StUtzstelle t n _i 

— B(x n ) ein zweiter Wert an der StUtzstelle x n , 
5 bezeichnet wird, und 

— bei dem durch iterative Losung des Naherungsprozesses die 
Storung Xx n + i ermittelt wird. 

2. Verfahren nach Anspruch 1, 

10 bei dem die Storung Rauschen ist, dem das System unterliegt. 

3. Verfahren nach Anspruch 1 oder 2, 
bei dem die Storung rein additiv ist. 

• 6 

> . I * • * 

. 4 f 

15 4. Verfahren nacfi einem der AnsprUche 1 bis 3, 
bei dem die Storung pfadweise ermittelt wird. 

5. Verfahren nach einem der Anspruche 1 bis 4, 

bei dem ein-stetiger Naherungsprozeii {x s ;s e [0, t]} ermittelt 

20 wird, indem ciie ermittelten Realisierungen des Naherungspro- 
zesses interpoliert werden. 

6. Verfahren nach Anspruch, .4, und 5, 

bei dem die ermittelten Pfadd mit einem stiatistischen Verfah- 
25 ren analysiert werden. 

7. Verfahren nach einem der AnsprUche l,bis 6, 
bei dem das System eine elektrische "Schaltung ist. -"' 

30 8, Verfahren nach einem der Anspruche 1 bis 6, 

bei dem das System ein mechanisches Mehrkorpersystem ist. 

9. Verfahren nach einem der Anspruche 1 bis 6, 
bei dem das System ein physikalisches System ist. 

35 • •• 

10. Verfahren nach einem der AnsprUche 1 bis' r 6, 
bei dem das System ein- "clremisches- System ist. 
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11. Verfahren nach einem der Anspruche 1 bis 6, 

bei dem das System ein physikalisch-chemisches System ist. 

5 12. Vorrichtung zur Ermittlung einer Storung eines -techni- 
schen Systems, welches einer Storung unterliegt, 
bei der eine Prozessoreinheit vorgesehen ist, die derart ein- 
gerichtet ist, dafi 

— das System mittels eines impliziten stochastischen Diffe- 
10 rentialgleichungssystems beschrieben wird, 

— eine naherungsweise Losung des Systems ermittelt wird, in- 
dem ein diskreter Naherungsprozefi realisiert wird, 

— der diskrete Naherungsprozeil gemafl folgender Vorschrift 
realisiert wird: 

15 

— (c + ha 2 c) • X Xn + 1 - {- (l - y)c + h(l + yai - a 2 )c} • X Tn + 

+{- yc + h(y - aiy)c} • X Xn ^ 1 + 
+ha 2 s(x n + 1 ) + h(l + yai - a 2 )s(x n ) + 
+h(y - aiyjsfxn.x) + 
+B(x n ) • AWn + yB(x n ) • AW n _! 

wobei mit 

— C eine erste Matrix, 

— ai, <x 2 , y vorgegebene Parameter aus dem Intervall [0, 1], 

T 

20 — h: = — , eine Schrittweite in emem Ausgangszeitintervall 
N 

[0, T], wobei T ein vorgegebener Wert ist, welches in N 
Teilintervalle unterteilt ist, 

— G eine zweite Matrix, 

— X T ^ eine Realisierung des Naherungsprozesses an einer 

l n + l 

25 StUtzstelle x n + 1 , 

— X Xn eine Realisierung des Naherungsprozesses an einer 

StUtzstelle x n , 

— X T ^ eine Realisierung des Naherungsprozesses an einer 

c n — 1 

StUtzstelle x n _i, 
30 — s(x n + 1 ) ein erster Wert an der StUtzstelle x n+1 , 
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— s(x n ) ein erster Wert an der Sttttzstelle x n , 

— s(x n _i) ein erster Wert an der Stiitzstelle x n _i, 

— ein Differenzwert A%_i: = W Tn - w x n _ x zwischen einem zwei- 
ten Wert W Xn an der Sttttzstelle x n und dem zweiten Wert 

5 w x n _i an der Sttttzstelle x n _i 

— B(x n ) ein zweiter Wert an der Sttttzstelle x n , 
bezeichnet wird, und 

— bei dem durch iterative Losung des Naherungsprozesses die 
Storung X Xn + 1 ermittelt wird. 

10 

13. Vorrichtung nach Anspruch 12, 

bei der die Storung Rauschen ist, dem das System unterliegt. 

14. Vorrichtung nach Anspruch 12 oder 13, 
15 bei der die Storung rein additiv ist. 

15. Vorrichtung nach einem der Anspruche 12 bis 14, 

bei der die Prozessoreinheit derart eingerichtet ist, dafi die 
Storung pfadweise ermittelt wird. 

20 

16. Vorrichtung nach einem der Ansprttche 12 bis 15, 

bei der die Prozessoreinheit derart eingerichtet ist, dafi ein 
stetiger NaherungsprozeB {x s ; s <= [0, t]} ermittelt wird, indem 

die ermittelten Realisierungen des Naherungsprozesses inter- 
25 poliert werden. 

17. Vorrichtung nach Anspruch 15 und 16, 

bei der die Prozessoreinheit derart eingerichtet ist, dafi die 
ermittelten Pfade mit einem statistischen Verfahren analy- 
30 siert werden. 



18. Vorrichtung nach einem der Anspruche 12 bis 17, 

bei dem das System eine elektrische Schaltung ist. 

35 19. Vorrichtung nach einem der Ansprttche 12 bis 17, 

bei dem das System ein mechanisches Mehrkorpersystem ist. 



WO 99/48030 



PCT/DE99/00239 



37 

20. Vorrichtung nach einem der Anspruche 12 bis 17, 
bei dem das System ein physikalisches System ist. 

21. Vorrichtung nach einera der AnsprUche 12 bis 17, 
bei dem das System ein chemisches System ist. 

22. Vorrichtung nach einem der Anspruche 12 bis 17, 

bei dem das System ein physikalisch-chemisches System ist. 
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FIG 1 



Speicherung der Schaltung in Rechner in 
Form einer Schaltungsbeschreibungssprache 



modifizierte Knotenanalyse der Schaltung 



•101 
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Ermittlung der Storung Xx n +1 durch Losung des 
Naherungsprozesses iterativ fur Zeitpunkte tp+i 
n =0.1 N-1 

-{C+h a2 G}.X T n+1={l 1 -Y) c+n ( 1+ W- a 2) G > X Tn 
+ {-yC+h(y-a-|Y) G}-Xt n _i 

+ha2$ (-c n +l)+h(1+yai-a2) s (xn)+h (y-a-iy) s (TrM) 
+B(xn)-AW n+Y B(Tn)-AWn-1 

103 
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